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Abstract 

The constraints on the minimal supergravity model from the b ^ s'y decay are stud- 
ied. A large domain in the parameter space for the model satisfies the CLEO bound, 
BR(6 —* 57) < 5.4 X 10~^. However, the allowed domain is expected to diminish signifi- 
cantly with an improved bound on this decay. The dependence of the 6 — * 57 branching 
ratio on various parameters is studied in detail. It is found that, for At < and the top 
quark mass within the vicinity of the center of the CDF value, m^°^^ = 174±17 GeV, 
there exists only a small allowed domain because the light stop is tachyonic for most of 
the parameter space. A similar phenomenon exists for a lighter top and At negative when 
the GUT coupling constant is slightly reduced. For At > 0, however, the branching ratio 
is much less sensitive to small changes in mt, and ao- 
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The extensive analyses of the high precision LEP data in the last few years have 
indicated that the idea of grand unification is only valid when combined with supersym- 
metry [1]. One of most promising and most studied models is the minimal supergravity 
model (MSGM) [2]. Supersymmetry (SUSY) is naturally and softly broken by a hidden 
sector. In addition to Yukawa couplings, the gauge coupling constant, ckg, and the unifi- 
cation scale, Mqut — lO"*^^ GeV, there are only five free parameters in this model: the four 
soft breaking terms [the universal scalar mass, mo, the universal gaugino mass, mi/2, the 
cubic scalar coupling, Aq, and the quadratic scalar coupling, Sq], and a supersymmetric 
Higgs mixing, iiq. The common approach to constraining the GUT model is to utilize the 
renormalization group equations (RGE) to make contact with physics at the electroweak 
scale. Mew [3]- Remarkably, the evolution of RGEs from Mqut to Mew produces a 
Higgs potential with a negative m|^^ if the top quark is heavy, signaling a spontaneous 
breaking of the electroweak gauge symmetry. As a consequence of this radiative breaking, 
two constraints arise which relate the GUT parameters to the electroweak parameters. 
One may then eliminate two GUT parameters, Bq, and in favor of two electroweak 
parameters: the Higgs VEV ratio tan/? = v^jvx^ and the Z boson mass, Mz- Therefore, 
the low-energy physics depends only on four parameters 

mo, m~g, At, tan/?, (1) 

and the sign of fj, [since the renormalization group equations determine only //^] . Here, we 
have replaced mi/2 by the gluino mass, rrig, and Aq by At, their values at the electroweak 
scale. Mew- In most analyses, only the masses of the third generation of leptons and 
quarks are important. For small tan/3, one needs to retain only the top quark mass. Hence, 
the four free parameters in eq. (1), plus the top quark mass, mt, suffice to parametrize the 
MSGM. 

The constraints on the parameter space of the MSGM may be classified into three 
major categories. First, various theoretical considerations put stringent constraints on the 
parameter space. For example, the color SU(3)(^ group should remain unbroken when 
discussing the radiative breaking; fi^ should also stay positive to guarantee this breaking; 
all scalar particles must be non-tachyonic; and the allowed parameter space should be such 
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that theory remains in the perturbative domain. Some of these issues will be discussed 
in a separate paper [4] . Second, cosmological considerations and the proton stability also 
strongly constrain the model [3]. Thus, for SU(5)-type models, proton stability requires 
that tan/3 should not be too large, i.e., tan/?<10 [5]. There still exists a large domain in 
the parameter space which satisfies both the proton decay and the relic density bounds [6] . 
Third, there exist a vast amount of data from the electroweak physics in the low-energy 
domain. It is thus very interesting to use these data to constrain the MSGM. One of the 
interesting processes is the b ^ S'j decay. This decay is very sensitive to the structure of 
fundamental interactions at Mew, because its rate is of order G^a, while most other FCNC 
processes are of order G^a^. We shall study the constraint coming from the h ^ decay 
in the MSGM in this paper. The combined constraints, from h — > s7,proton decay,and 
relic density will be discussed elsewhere [7] . 

The recent CLEO II experiment gives the following measurement for the branching 
ratio of the exclusive B K(892)*7 decay [8], BR{B K(892)*7) = (4.5 ± 1.5 ± 
0.9) X 10~^, and an upper bound on the inclusive B Xg^ decay, BR(S — > Xg^) < 
5.4 X 10~^. Much work has been devoted to the determination of the ratio BR(S 
JC(892)*7)/BR(i? — ^ X.s'y) from QCD calculations [9-11]. An accurate calculation is of 
interest since the inclusive B decay is not as easy to study experimentally as the exclusive 
B decay. The earlier calculations [9] made use of the constituent-quark model and the 
2-point QCD sum rules, and the results for this ratio ranged from 0.05 to 0.40. More 
recently, there have appeared calculations based on some new developments in QCD. The 
method based on the 3-point QCD sum rules predicts a value of 0.17 ± 0.05 [10]. The 
method that incorporates the chiral symmetry for the light degrees of freedom and the 
heavy quark spin-flavor symmetries for the heavy quarks gives a value of 0.09 ± 0.03 for 
this ratio [11]. Although these results are in accord with the CLEO II measurement, 
more study is needed to fully utilize the CLEO result. Because the b quark mass is much 
greater than the QCD scale, the dominant contributions to the inclusive B decay are from 
the short-range interactions. Barring a large interference between the long- and short- 
range contributions, one may model the inclusive B transitions, such as B ^ Xg'y and 
B — > XceVf,, by the decay of a free b quark to a free light quark, such as a free s or c quark. 
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In fact, a calculation based on the heavy quark effective theory shows that BR(S — > Xg^) 
agrees with the free quark results, BR(6 57) up to corrections of order 1/m^ [12]. Thus, 
the bound on BR(i? ^s7) transfers to an upper bound, BR(6 ^ 57) < 5.4 x IQ-^. It 
is a common practice to use the ratio defined as 

BR(6^g7) _ BR(^^X,7) 

to constrain various models, utilizing the well-determined value of (10.7±0.5)% for BR(S 
XceVg). The advantage of using i?, instead of BR(6 — > 57), is that the latter is dependent 
upon and certain elements of the Kobayashi-Maskawa matrix, while the former only 
depends on 2; = mc/wb, the ratio between the c and h quark masses, which is much better 
determined than both masses, i.e., z — 0.316 ± 0.013 [13]. 

The ratio R defined in eq. (2) has been calculated as [14] 

/(.)[l-AM^/(,)] 

where r] = as{Mz) / asirrih) = 0.548. Here, I{z) = 1 - 8z'^ + 8z^ - z^ - 2Az'^\nz is 
the phase-space factor, and f{z) = 2.41, a QCD correction factor, for the semileptonic 
process, b —>■ ceve- Ci^Mz) and Cs{Mz) are the coefficients of the photonic and gluonic 
penguin operators for the hs transition at the electroweak scale. These coefficients are 
model dependent and sensitive to the underlying fundamental interactions at Mz- For the 
standard model (SM), only the penguin diagram induced by the W-t loop contributes to 
Ci{Mz) and Cs{Mz), whereas, for the MSGM, many susy particles contribute. C2{Mz) is 
a coefficient coming from a mixing between the photonic penguin operator and many four- 
quark operators present at Mz- The form and number of these four-quark operators differ 
depending on the model. Fortunately, for the MSGM, they are the same as those for the 
SM [15]. The calculation of C2{Mz) is an involved procedure and a number of evaluations 
exist [14,16,17]. We shall use the results of Ref. [17], which takes into account of the full- 
leading-order logarithmic contributions, i.e., C2{Mz) — Ylii=i'^i''l — —0.1795, where ai 
satisfy = [since at the electroweak scale, only the photonic penguin diagram has a 

contribution to the 6 ^ 57 decay] . The numerical values for ai and hi have been a matter 



of debate, in the sense that different calculations give different results. Nevertheless, the 
over all effect on C2{Mz) is insignificant, the difference being only about 1%. A comment 
on the accuracy of eq. (3), however, is in order. There exist several uncertainties in using 
this equation to calculate BR(6 — > sj). For instance, eq. (2) is based on the spectator 
quark model. The error in determining the strong interaction constant, as{Mz), is still 
large. The most significant uncertainties come from the absence of a complete evaluation of 
the next-to-leading short-distance QCD corrections to BR(6 — > 57), causing about a 25% 
inaccuracy. To better determine the theoretical predictions for BR(i? Xs'j), it would 
be necessary to calculate certain three-loop mixings and two-loop penguin diagrams [18]. 

In the rest of this paper, we will concentrate on the MSGM predictions for BR(6 
57). As mentioned above, this model contains many particles not present in the standard 
model. Thus, besides the W boson contributions, there exist the penguin diagrams induced 
by the charged Higgs bosons, the chargino-squarks, the gluino and the neutralinos. The 
coefficients C7{Mz) and Cs{Mz) were calculated for the minimal supersymmetric model 
(MSSM) in Ref. [19]. (We have rederived these coefficients confirming their results.) It is 
found that the gluino and neutralino contributions are small compared to other sources. 
We will hence ignore their contributions below. A large contribution from the charged 
Higgs boson in the MSSM was found in Ref. [20] and these authors thus concluded that a 
slight improvement in the experimental bound on BR(6 S7) will exclude the search for 
the charged Higgs boson via the t — > bH'^ decay channel. However, these papers did not 
consider the chargino-squark penguin diagrams. As it turned out, although the charged 
Higgs boson enhances the standard model amplitude, the chargino-squark loops may con- 
tribute to the amplitude constructively or destructively, depending on the parameters 
chosen. In fact, as shown in Ref. [21], in the exact supersymmetric limit, the coefficients 
for 6s7 and bsg transition operators, C'y{Mz) and Cs{Mz)., vanish exactly. Other papers 
on the 6 — > 57 decay in SUSY models can be found in [22]. We will follow the notations 
of Ref. [21], and assume that the first two generations of squarks arc degenerate in mass. 
We then expect the contributions from these degenerate squarks to BR(6 — > 57) to be 
small, since their masses are proportional to mo^l TeV. On the other hand, the scalar 
top squarks, ti and t2, are badly split in their masses, due to the large top mass, implied 
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by the recent CDF data [23] . The stop mass matrix is given by 



m? nit {At + n cot (3) 

2 



nit {At + iJ, cot l3) mj 



(4) 



where and m?^ are given in Ref. [24]. Thus, the hght stop mass is, 

1 = 1 (""I + ™L - ^("^i - mlj^ + 4mf{At + /icot /3)2) . (5) 



ti 



One can demonstrate that, for a large portion of the parameter space, the light stop 



(mass)^ may turn negative, signalling either the breaking of the color SU(3)^ group or the 



existence of tachyons in the theory. The requirement that m| be positive is very stringent 
and could eliminate a considerably large domain in the parameter space [4]. 

Our strategy is to use the one-loop RGEs to calculate the mass spectrum of the 
model relevant to 6 — > 57 [all the up-type squarks, the charginos, and the charged Higgs 
bosons]. We then use this spectrum to evaluate the ratio R via eq. (3) and eq. (2) to obtain 
BK{B — > Xg^) — R ■ BR(i? Xf-eVg), and compare the results with the CLEO II bound. 
Of the five parameters, rrit is restricted by the CDF bound, mf"^'^ ~ 174 ± 17 GeV [23]. In 
our analysis, mt is the running top mass at M^, which is related to the pole mass by [25] 

mf'^ = ruMl + ^-^^ + ll{^^n (6) 

where p* ~ ht{mt) / ht{M z) is the ratio of the top Yukawa couplings at mt and Mz [We 
assume that the Higgs VEVs, vi and V2, and hence tan/3, do not change significantly 
between mt and Mz\- The pole mass given by eq. (6) is about 3 — 5% larger than the 
running mass at Mz- The naturalness condition restricts mo and mg to be less than 
0(1 TeV), and we allow them to lie between 100 GeV and 2 TeV. We parametrize At in 
units of mo, and restrict |^t/mo| < 2.0. Notice also that the top mass is very close to 
its Landau pole, mf""^'*" = Csin/3, where C ~ 195 GeV. This implies that the error m 
as{Mz) plays an important role in the RGE analysis of the spectrum. Here, instead of 
varying as{Mz)-i we let ac vary between 1/24.11 and 1/24.5 with a fixed Mqut — 10^^- 
since a slight change in Mqut does not affect the coupling unification as significantly as a 
change in aa- It turns out that aa = 1/24.11 gives, at the one- loop level and for m^ = M^, 



the best fit to as{Mz) = 0.118, a2{Mz) = 0.03358 ± 0.00011, and ai{Mz) = {5/3)aY = 
0.016985 ±0.000020 for the Mqut cited above, and the two-foop correctfons are small [4]. 
Decreasing ao corresponds to decreasing as{Mz). We will discuss below the consequences 
of varying ao- 

We have surveyed a large domain in the parameter space described by eq(l). The 
branching ratio depends importantly on the values of the parameters in this 4-dimensional 
space. Significant deviations from the SM value can occur in certain regions of the pa- 
rameter space. Characteristically the region where large deviations from SM arise occur 
when mo, rrig are much smaller than their naturalness limits and tan/3 gets large, i.e., 
typically larger than 10. Specifically in this region of the parameter space BR(6 — > 57) 
can be significantly below its SM predictions, and for certain points, an almost perfect 
cancellation is observed. This is because the chargino-squark penguin diagrams contribute 
destructively to the total amplitude at these points, with a coefficient ~ 1/ cos /3 ~ tan (3. 
This destructive interference between various sources in MSSM has also been observed 
previously [26] [the symmetric distribution of the branching ratios around the SM values 
found in Ref. [26] is because those authors allowed tan P to be as large as 60] . Similar de- 
viations from the SM value can also occur for smaller tan/3, although this is less frequent. 
For tan/3 < 10.0, the deviations of BR(6 — ^ S7) from the SM values are less dramatic. 
In this region, the current CLEO bound is not stringent enough to strongly constrain the 
MSGM. However, with a moderate, e.g., about 30%, improvement in the CLEO bound 
constraints on the model will emerge. 

Figures 1, 2, and 3 show plots of BR(6 S7) as a function of the light chargino 
masses, w^^^, and the soft supersymmetry breaking parameter, mo, for tan/3 = 5.0, 
lAt/niol = 0.5, mt = 160, 170 GeV, and a^^ = 24.11, 24.5. In these graphs, all the 
masses are in units of GeV. The graphs for rrit = 150 GeV are similar to Figures 1, 2, 
and 3, except that the branching ratios are smaller. We also impose a phenomenological 
lower bound on the light chargino mass, i.e., rriy^^^ > 45 GeV. These figures are character- 
ized by four parameters, mt, ct^^, and the signs of At and /i. The figures labeled by 'a' 
have At < 0, /J, > 0, those by 'b' have A^ < 0, /j, < 0, those by 'c' have At > 0, ij, > 
and those by 'd' have At > 0, /j, < 0. These figures contain the following results. (1) 
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The b ^ branching ratio at the points at which At and have the same sign is in 
general larger than that at the points where At and fi have the opposite sign. This can 
be seen by comparing Figures la and Id with Figures lb and Ic. This is a generic feature 
of the MSGM, because eq. (5) gives a smaller light stop mass when At and n have the 
same sign. The gaps in the lines for mo = 1000 GeV in Figure la and for mo = 600 and 
1000 GeV in Figure lb are due to the light stop turning tachyonic. These gaps occurs 
only for At negative for the region in parameter space we have studied. Generally, the 
constraint that the light stop remains non-tachyonic is most stringent for At < 0. (2) The 
6 — s> S7 branching ratio increases when rrit increases. The dependence on nit is correlated 
with other parameters. The parameters for Figure 1 (m^ = 160, mf°'^ ~ 165) and those 
for Figures 3 (m^ = 170, m^"^" ~ 175) differ only in rrit. Remarkable is that a slight change 
in mt (about 6%), for At negative (Figures la and lb vs. Figures 3a and 3b) eliminates 
a large part of the allowed parameter space. For example, while the values of mo ranging 
from 100 to 1000 GeV are all allowed for rrit — 160 GeV, for both positive (la) and nega- 
tive (lb) fj,, the only allowed values of mo are below 200 GeV for rrit = 170 GeV for both 
positive (3a) and negative (3b) However, if At is positive (Figures Ic, Id, 3c, and 3d), 
the same change in rrit does not significantly affect the allowed parameter domain. The 
reduction in parameter space can again be explained by the light stop turning tachyonic. 
(3) Most interesting is the a.c dependence of the branching ratio. To see this effect, let us 
compare Figures 1 [cu^^ = 24.11, for which as{Mz) = 0.118] with Figures 2 [a^^ = 24.5, 
for which as{Mz) = 0.113]. For At negative, we again find that a larger part of the pa- 
rameter space are excluded (la vs. 2a, and lb vs. 2b), while for At positive, the allowed 
parameter space remains almost the same (Ic vs. 2c, and Id vs. 2d). Although the light 
stop turning tachyonic is the reason for this, the physics involved is quite different from 
the above. Qualitatively, one can attribute this phenomenon to the fact that ht is very 
close to its Landau pole. Namely, a reduction in ao modifies various form factors defined 
in Ref. [24], making ht closer to its Landau pole. This in turn is reflected in the light 
stop turning tachyonic. The combined effects of simultaneous change in mt and aa is 
very dramatic — the only allowed parameter space is mo = 100 GeV for mt = 170 and 
— 24.5 when At is negative! This is because the simultaneous changes in both mt 
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and ao add up almost multiplicatively to aggravate the closeness to the Landau pole [4] . 
For At positive, the change is not as large, and the allowed domain is still large. Similar 
reduction of the allowed domain also exists for At with other negative values - the more 
negative At is, the smaller the allowed domain remains. 

In conclusion, we have performed a detailed study of the constraints from the h ^ 
decay on the MSGM. There are regions of the parameter space where the branching ratio 
exceeds the CLEO II bound, and this region is excluded . However, there still exists a large 
domain that satisfies the CLEO bound. A more accurate determination of the branching 
ratio would further constrain this model. An interesting result of the analysis is that very 
little allowed domain of the parameter space was found for At < Q and mt in vicinity of 
the CDF central value of 174 GeV for < 0. Thus the allowed domain resides mostly in 
> region. 

This work was supported in part under NSF Grant Nos. PHY-916593 and PHY- 
9306906. 
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Figure Captions 

Fig. 1 The branching ratio, BR(6 57), plotted with respect to the hght chargino mass, 
m^y^ > 45 GeV, for tan/3 = 5.0, the running mass nit = 160 GeV, |At/mo| = 0.5, 
and a^^ = 24.11, and for (a) At < 0, /i > 0; (b) At < 0, // < 0; (c) At > 0, > 0; 
and (d) At > 0, /j, < 0. All masses are in units of GeV. The standard model gives a 
branching ratio, BR(6 — > S7) = 3.55 x 10~^. The discontinuity in the mo = 1000 GeV 
line for (a) and in the mo = 600, 800 and 1000 GeV lines for (b) is due to the light 
stop turning tachyonic. 

Fig. 2 Same as Fig. 1 for aa = 1/24.5. The discontinuity in the mo = 1000 GeV line for (a) 
and for (b) is also due to the light stop turning tachyonic, as in Figure 1. 

Fig. 3 Same as Fig. 1 for the running mass rrit = 170 GeV. The standard model gives a 
branching ratio, BR{b 57) = 3.66 x 10"^. 
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